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Using a general model we show that under a common external forcing, the species with a quadratic 
saturation term in the population dynamics first undergoes spatial synchronization and then extinc- 
tion, thereby avoiding the rescue effect. This is because the saturation term reduces the synchro- 
nization time scale but not the extinction time scale. The effect can be observed even when the 
external forcing acts only on some locations provided there is a synchronizing term in the dynamics. 
Absence of the quadratic saturation term can help the species to avoid extinction. 
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Consider two important phenomena concerning popu- 
lations of different species. First is the spatial synchro- 
nization of populations of a species. Many examples of 
spatially synchronized populations have been observed in 
nature ^Q^^A-,^3eIL^2x jjj 0> 113 • These include 
synchronization of vole populations by predatory birds 
W , synchronization of caribou and musk oxen by climate 
[i 12j and synchronization of lynx populations probably 
by climate [g and/or dispersalja. Several other exam- 
ples are documented in Ref. [lOj. The second impor- 
tant phenomena is the extinction of species. More than 
99% of the species that ever existed on the surface of the 
earth are now extinct. One example of such extinctions 
is the statistically homogeneous K — T extinction inten- 
sities observed for marine molluscs on a global scale |0| . 
(In this letter, extinction will mean extinction on the 
global scale.) Such global extinctions are still a puzzle. 
It is possible that when a species is under threat it may 
survive in some isolated locations and afterwards lead 
to the revival of the species. This is the 'rescue effect' 
[T^. Various factors affecting extinction such as migra- 
tioUj chaos, noise etc have been discussed in the literature 
[iSl llq . It is believed that spatial synchronization can 
promote global population extinctions 0,0. However, 
there is no clear understanding of the relation between 
spatial synchronization and extinction and whether they 
will always co-exist. This underscores the need for a gen- 
eral theory of spatial synchronization and extinction of 
populations under external forcing which will clarify this 
relation. 

In this letter we investigate the time scales of spatial 
synchronization and extinction. Using a general model 
of population dynamics under a common external forcing 
we show that the saturation term (decay term) decides 
these time scales. If the dynamics has a quadratic satu- 
ration term, the species first undergoes spatial synchro- 



nization and then extinction. Here, extinction takes place 
almost simultaneously and synchronously in all the loca- 
tions. Thus, the rescue effect can not take place. On 
the other had, the species which do not possess the spe- 
cific saturation term in the dynamics, will show a natural 
resistance towards extinction through the rescue effect. 

We present the details of our argument by consider- 
ing a general model of population dynamics. We follow 
the experimental set up of Ref. 1] that considered 28 
enclosed vole populations that were fenced to prevent 
predatory mammals and vole dispersal. It was found 
that the vole populations synchronized due to predatory 
birds. Following this experimental set up, let Pi (t) denote 
the population of a species at i—th patch, i = 1, ... TV, at 
time t and let Q{t) denote an external variable (e.g. me- 
teorite impacts, volcanic eruptions, predator population, 
climate etc.) which interacts with the population of the 
species at different patches. The coupled dynamics can 
be written as 

^ = h{P,{t)) + emim),Q{t))+I, (la) 

^ = f2m)) + ^T.92im),Qit)), (lb) 

i 

where /i and /2 represent the uncoupled dynamics, gi 
and (72 represent the interactions and Ci and £2 are inter- 
action constants. The term I represents the interaction 
terms between populations of species in different patches 
and is not important for our basic argument. We will 
consider its effect later. 

We are interested in the spatially synchronized state, 
Pi{t) — . . . — P/v(i) = P{t)- Linear stability of this state 
can be analysed using the following Jacobian J of the 
vector field defined by P^'s and Q. 
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The eigenvectors of J split into two orthogonal sub- 
spaces A and B 0. The subspace A has dimension 
two and it defines the synchronization manifold. The 
subspace B has dimension — 1 and it defines the trans- 
verse manifold. The eigenvectors of B arc of the type 
ctt — (cki, . . . , QfAT, 0)-'", Eflj^Qfi = and the eigenvalues 
are N ~ \ fold degenerate and are given by -§p{fi + tigi)- 
We note that under time evolution the two subspaces A 
and B do not mix with each other and the subspace B has 
the same eigenvectors for all the time. Thus, it is possi- 
ble to take the time average of the degenerate eigenvalue 
of B to obtain the transverse Lyapunov exponent. For 
the stability of the spatially synchronized state we re- 
quire the transverse Lyapunov exponent to be negative. 
Imposing this, we obtain the condition for the stability 
of the synchronized state as 

(^(/i(^) + ei5i(^,O)))<0, (2) 

where ( ) represents the time average. 

We now return to the problem of synchronization and 
extinction. Near extinction. We analyse the problem by 
retaining the lowest order terms in the Taylor series ex- 
pansion of various functions in Eq. Hla|l in terms of the 
population Pi{t). 

h{P,) = aP,-hP^ + 0{Pf). (3) 

The first term in the expansion is a growth term and 
a > 0. The second term is a saturation term if 6 > 0. 
Neglecting the higher order terms, we get the stable so- 
lution Pi = h/ a. If 6 < 0, then we must include higher 
order terms in Eq. Q to get a stable solution. The in- 
teraction function gi to lowest order in Pi can be written 
as —Pih{Q) where h{Q) is some function of Q. 

Let us first consider the case 6 > 0. The condition Q 
for the stability of synchronized state now becomes 

K = {a-hP -eih{Q)) <Q. (4) 

For extinction the forcing must be able to compensate 
the growth and the condition for extinction is 

\e = {a-eih{Q)) <Q. (5) 

This condition can also be obtained by considering the 
stability of P = state. Comparing the conditions Q 
and (O, we find that as (a — eih{Q)) starts decreas- 
ing, the condition Q will be satisfied before the extinc- 
tion condition (O is satisfied. If both the synchroniza- 
tion and extinction conditions are satisfied then the time 
scale associated with synchronization (r^ = 1/ |As|) will 
be less than the time scale associated with extinction 
(re = l/|Ae|). Thus, we conclude that the populations 
in different locations will synchronize before the extinc- 
tion of the species. They will remain synchronized as the 
populations at different patches start decreasing. Hence 
the extinction of populations in different patches will take 
place almost simultaneously. 



We demonstrate that spatial synchronization precedes 
extinction using a simple prey-predator model |l8j |. 
For this model, the different functions in Eq. (I) are 
given by h{P) = aP - bP\g^{P,Q) = -g2{P,Q) = 
~PQ,h{Q) = -u{Q - Q*). We allow the predator to 
maintain a low equilibrium level Q = Q* even when its 
usual prey, P, is rare ^|. For the above model, the 
synchronization condition becomes (a — 2bP — tiQ) < 
and the extinction condition is given by (a — eiQ) < 0. 
In Figure 1(a), we plot populations of different patches 
as a function of time starting from random initial popu- 
lations. The parameters used are: a — 0.5, b — 50.0, u — 
0.1, ei = 4.8, €2 = 1.0, N = 100, Q* = 0.5. We see that 
the populations of different patches synchronize and then 
are driven to extinction. To better understand the time 
scales involved, we plot the following two parameters as a 
function of time in Figure 1(b): Synchronization param- 

eter S = {2/N{N - 1)) Y.f=i Y.]~=\{P^ " P3Y which mea- 
sures the mean square deviation between pairs of pop- 
ulations and extinction parameter E — {^/N)J2iLiPi 
which measures the mean square populations. We ob- 
serve that initially the synchronization parameter S 
(solid line) goes to zero with a rate greater than that 
of the extinction parameter E (dotted line) . As the pop- 
ulations become very small, and Ag [Eqs. Q and jSJ] 
become nearly identical and the rates of decrease of S 
and E become nearly equal as can be seen from Figure 
1(b). 

In the dynamics of the populations we have neglected 
the effect of intra-species interactions or diffusion within 
the populations in different patches (/ in Eq. (|la|l ). It 
is easy to sec that these effects do not affect the con- 
clusion of spatial synchronization before extinction. For 
example, consider Eq. Hla() with 

N 1 ^ 

where the first term gives a mean field type interaction 
while the second term represents the interaction between 
different populations, e.g. /i2 = PiPj- The conditions for 
the stability of the synchronized state and the condition 
for extinction now become, respectively, 

([F(P,X,g)]^^p) <0, (7a) 
{[F{P,X,Q)]^^P^,)<0, (7b) 

where F{P, X, Q) = ^[/i(P) -f h2{P, X) - h^iX, P) + 
eigi{P, Q)]. Note that hi does not contribute to the above 
conditions and /i2 does not contribute if it is symmetric 
in its arguments, e.g. h2{X,Y) = h2{Y,X). We see 
that our conclusion about spatial synchronization before 
extinction is still valid. Other effects like migration or 
diffusion can be treated in a similar fashion. Again our 
basic argument remains valid. 

Our basic finding of synchronization before extinction 
can be tested in an ecological experiment similar to those 
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described in Ref. [jj] . The food source of the vole popu- 
lation can be decreased progressively to see whether the 
synchronization persists and whether the extinction is al- 
most simultaneous. Similar experiments could be carried 
out with other populations such as insects etc. 

In the above argument we have neglected the effect of 
higher order terms in the Taylor series expansion of vari- 
ous functions appearing in Eqs. (|lal) . When higher order 
terms are retained two possibilities arise. First chaotic 
attractors can occur and secondly there can be multiple 
stable solutions . It is believed that when the isolated 
patches are individually chaotic with a weak coupling 
then it can lead to asynchrony between different patches 
thus preventing global extinction [l^ . However, we find 
that such chaotic solutions do not prevent spatial syn- 
chronization when the effect of common external forcing 
is included. It is well known in the nonlinear dynamics 
literature and also as shown above, that a common forc- 
ing, including common external noise, can synchronize 
chaotic systems [23, 0| . 

When there are multiple stable solutions, it is possi- 
ble that depending on the initial conditions populations 
in different patches may converge to different stable so- 
lutions. In this case extinction is difficult due to rescue 
effect and neither is there a spatial synchronization. 

So far we have neglected the effect of variation in pa- 
rameters from patch to patch and also the effect of lo- 
cal noise. Both these factors, if they are large, can lead 
to spatial asynchrony. However, we have verified that 
small parameter variations and noise do not affect our 
conclusions. The effect of small parameter variations is 
demonstrated afterwards in Fig. 2. 

We now discuss the case 6 < in Eq. Q . It is easy to 
see that in this case the time scale for extinction will be 
less that that for synchronization. Clearly, spatial syn- 
chronization cannot take place before extinction. Thus, 
it is is possible that the rescue effect can prevent extinc- 
tion. In Fig. 1(c), we show the poputions of different 
patches as a function of time for b = —4.0 We see that 
the patches do not show any spatial synchronization and 
the species may survive due to rescue effect if the ex- 
ternal forcing is switched off after some time. Fig. 1(d) 
shows the synchronization and extinction parameters, S 
and E, as a function of time for b = —4.0. We see that 
initially the rate of decrease of E is greater that that of 
S and afterwards the two rates become almost equal. 

Thus, we see that the parameter b can be treated as a 
measure of the resistance of a species towards extinction. 
Smaller the value of b, more is the resistance. The second 
term in Eq. (|3Jl corresponding to parameter b represents 
interaction between two members of a species. The pa- 
rameter b is in general positive due to competition be- 
tween members. This is also reflected in various popu- 
lation models used in the literature which are known to 
give good fit for experimental observations 0, [S^, ■ 
However, a high degree of cooperation between the mem- 
bers may be able to make b negative and the species more 
resistant to extinction. 



We now consider a situation of great practical impor- 
tance. It is easy to establish that if the coupling param- 
eter ei has a small variation, it does not alter the above 
conclusions. However, it may happen that some patches 
escape the effect of the external forcing i.e. ei = for 
these patches. We now show that our conclusions based 
on the parameter b are still valid provided there is some 
synchronizing interaction among the patches. Let N2 
patches escape the effect of external forcing and the re- 
maining A^i — N — N2 patches be affected by the forc- 
ing. We choose the interaction term I — J2jiPj ^ Pi)- 
In this case, by using an argument similar to the one 
used to show the stability of the synchronized state for 
Eq. (la), it is possible to show that we get a two-cluster 
synchronized state. The A'^i patches synchronize to one 
value of the population say Pi and the remaining to an- 
other value say P2- The conditions for the two cluster 
synchronized state are — bPi — d — eih{Q)^ < and 
(a — 6P2 — d) < for the two clusters respectively. The 
difference AP — P2 — Pi evolves as 

^ = [a - 6(Pi + P2) - d\AP - eMQ)Pi (8) 

When the extinction condition for Ni cluster is satisfied 
i.e. (a — eich{Q) — d) < and 6 > 0, the synchronization 
of Pi will precede that of extinction due to the b term. 
The cluster P2 will also synchronize due to b term if P2 is 
large. Now as Pi becomes smaller the second term on the 
RHS of Eq. © becomes small and if [a-6(Pi + ^2) -rf] < 
then P2 will start decreasing in some sort of generalized 
synchrony with Pi . 

In Fig. 2 the time evolution of populations in differ- 
ent patches is shown when only 50% of the patches in- 
teract with the external forcing. We also introduce a 
5% patch to patch variation in all the parameters about 
their respective mean values. For b positive (Fig. 2(a)) 
we initially see the formation of the two-cluster synchro- 
nized state. Following this, the iVi cluster rapidly decays 
and is closely followed by the N2 cluster. For b negative 
(Fig. 2(b)) the populations again separate into two dis- 
tinct groups. The A^i cluster shows a rapid decay but 
the A'2 cluster shows a very small decay. This small de- 
cay comes from the pulling down effect of Eq. We 
note that as the A''! cluster becomes extinct the effect of 
Eq. © will also vanish and the N2 group will start its 
independent evolution thus escaping extinction. 

In this letter, we have established a clear connection 
between extinction and spatial synchronization of popu- 
lations. Under reasonably general conditions with exter- 
nal forcing we showed that spatial synchronization pre- 
cedes extinction when the parameter b > thus prevent- 
ing the rescue effect. On the other hand, for 6 < 0, the 
species can show a natural resistance to extinction. These 
conclusions are valid even if the external forcing acts only 
at some locations provided there is some synchronizing 
interaction between the populations. Clearly in mass ex- 
tinction events where there is a strong common external 
shock, the above conclusions should be valid. Even in 
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other situations, we expect our general conclusions to 
hold since they are based only on the parameter b of the 
local dynamics. 
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FIG. 1: This figure demonstrates the interplay between syn- 
chronization and extinction in a simple prey-predator model 
fl3 . (a) and (c) show the populations of different patches as 
a function of time for b positive and negative respectively, (b) 
and (d) show the synchronization and extinction parameters, 
S and E, (solid and dotted lines respectively) as a function 
of time. 



FIG. 2: The time evolution of populations in different patches 
is shown when only 50% of the patches interact with the ex- 
ternal forcing, (a) b = 50.0; (b) b = -4.0. 



